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^ : Abstract 
o 

Q ! A class of Stephani cosmological models as a prototype of non-homogeneous 

"^H I universe is considered. The non-homogeneity can lead to accelerated evolu- 

Qh' tion which is now observed from the SNIa data. Three samples of type la 

Q ■ supernovae obtained by Perlmutter at ai, Tonry et al. and Knop et al. are 

taken into account. Different statistical methods (best fits as well as max- 
! imum likelihood method) to obtain estimates of the model parameters are 

used. Stephani model is considered as an alternative to the concordance of 
IT/i ' ACDM model in the explanation of the present acceleration of the universe. 

j_j ■ The model explains the acceleration of the universe at the same level of ac- 

^ curacy as the ACDM model (x^ statistics are comparable). From the best fit 

analysis it follows that the Stephani model is characterized by higher value 
of density parameter Qmo than the ACDM model. It is also shown that the 
obtained results are consistent with location of CMB peaks. 

1 Introduction 

In the paper of Stelmach and Jakacka ^ it was suggested that the effect 
of acceleration of the universe can be driven by non-homogeneities in the 
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Stephani model [21 E] ■ We perform the statistical verification of this hipothe- 
sis using current available astronomical data. We use three different samples 
of supernovae where Perlmutter sample |lj is treated as the fiducial data set. 
For comparison both the Tonry et al. [HI (improved by Barris et al. [Hj) and 
Knop et al. [Ij are used. 

To find which model fits to the data in the best way we use the 
statistics. The parameters of the models are estimated using the maximum 
likelihood method. One-dimensional probability distribution function (pdf) 
over model parameters is presented to deepen statistical insight. 

Recently found accelerated expansion of our Universe jH] is explained in 
the literature by the presence of the A-term or the so-called quintessence 
matter satisfying the equation of state with negative pressure P . Observa- 
tion of type la supernovae (SN la) revealed that this form of energy called 
dark energy could dominate the present evolution of the universe. Therefore 
it is important to consider different candidates for dark energy. In the paper 
of Stelmach and Jakacka it was shown that non-homogeneity in spherically 
symmetric Stephani model can be treated as some kind of dark energy due to 
which the universe accelerates (we call this model "S-J model"). It is a sort 
of fictitious fluid for which density parameter f2non can be defined. Then dy- 
namics can be formally reduced to the FRW flat model with additional fluid 
satisfying the equation of state p = (— 2/3)p corresponding to the equation 
of state for topological defects in a form of domain walls fOJ El E] ■ 

The best fitting procedure applied to supernovae data as well as confi- 
dence level for the redshift-magnitude relation shows that Qnon should dy- 
namically dominate the present evolution of the Universe if we want to ex- 
plain observational data without cosmological term (see also |13j). 



2 Stephani Model in the S-J version 

Stephani model in the S-J version can be described by the generalized Fried- 
man equation [T] 

^' + ^ = ^' (1) 

where k{t) = f3R{t) is an effective curvature index which depends on time. 
(3 and a are constants {p = ap). In the above equation non-homogeneity is 
not explicitly present. This is due to the choice of the time parameter. 

Equation (1) is de facto a first integral of the Einstein equation which is 
a second order equation for the scale factor. It is obtained by assuming that 
the observer is placed at a symmetry center p{r ^ 0,t) = ap{t). We consider 
the universe in the neighbourhood of this center. 
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The relation 

kit)=(3Rit). (2) 
is a special form of a more generalized ansatz of a type 

= /3i?T(t),7 = const. (3) 

However, 7 = 1 is a simplest case, and our goal is to show that such exotic 
model can explain SN la data of Perlmutter even in this simple case. 

For further purposes it is convenient to write down Eq. (1) in a dimen- 
sionless form 



1 

■„2 



J2^.,ox-^'^^^'\ (4) 



X 

1=0 

where x = R/Ro is the so-called radius of the universe in units Rq. Here the 
differentiation is carried out over the dimensionless time parameter r = \HQ\t. 
The parameters f2jo = 0, 1) are defined 

-J -no 

where the subscript means that a quantity with this subscript corresponds 
to the present epoch 

pio = piix=l). (6) 
In our model = a and ai = —2/3 like for topological defects, hence 

- A -(-3/5) 

P'- = ^3^^ Pnon = — (7) 

where Ca is some constant. 

In the given parametrization our non-homogeneous Stephani model looks 
formally as a flat Friedman model with two noninteracting components. 



3 Stephani model as a hamiltonian system 

Because of interpretation of the S-J model in terms of FRW model with ad- 
ditional fictions fluid we can simply find hamiltonian formalism for dynamics 
of the system. A first integral of the Friedman equation can always be used 
for reduction of the motion of the system to a motion of a particle of unit 
mass in one dimensional potential (14j. If we define 

p,ff ^^p,(i?) = /(i?), (8) 
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then the potential is 



ViR)^-P^R' = ^R~--J— (9) 



and the Euler-Lagrange equation of motion reads 



dV 

and has a Newton-hke form. 

The Hamiltonian in the considered case has the form: 

n = ^ + v{R), (11) 

where we defined Pr = R 

Trajectories of the system he on a zero energy level Ti. = and a hamil- 
tonian constraint must coincide with the form of the first integral of the 
equation of motion (10). 

The advantage of having the hamiltonian function given by Eq (11) is 
that one is then able to make easily full classification of solutions in the con- 
figuration space P] . Moreover dynamics can be presented in two dimensional 
phase plane. 

The corresponding equations of motion describe a simple dynamical sys- 
tem consisting of a particle moving in one dimensional potential: 

x = y, 



where x = R, y = R and V{x) is given by (9), and j/^/2 -|- V{x) = is 
its first integral. Qualitative analysis of differential equation implies a shift 
from finding and analysing individual solution to investigating the space of 
all solutions. Certain properties (such as existence or absence of horizons, 
existence of singularities) are believed to be realistic if they can be attributed 
to a larger class of models within the space off all solutions (phase space). 

This approach offers a possibility of investigating the space of all possible 
solutions for the considered problem. Of course the system (12) possesses 
the first integral = which defines a family of algebraic curves (phase 
curves) on which the trajectories of the system lie. The phase portrait of the 
considered system for a) a = (dust) is shown in Fig. 1 and for b) a = 1/3 
(radiation) in Fig. 2. 
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The system is described by the equations x = y, y = —dV/ dx and the first 
integral is y'^/2 + V{x) = 0, where V{x) = (l/2)(/3a; — A^/x^^^") and we put 
(3 = —1. The phase domain is x > 0. The acceleration region is situated to 
the right of the saddle point, therefore for the Lemaitre-Eddington (loitering) 
universes the acceleration begins in the middle of a quasi-static stage. 

For comparison the FRW model with dust and cosmological constant 
is presented in Fig. 3. We can observe how non-homogeneity term can 
mimic cosmological constant term. In the particular cases we have: a) x = 
y^y = -{l/2){/3x + A^/x'^), h) X = y,y = -(l/2)(/5 + A^x^), where A = 
OA, P = —1, c) X = y,y = — (A/3)x — po/6x^, where A = 0.7, po = 0.3 with 
V{x) = -(A/3)x2 -po/6x. 

The phase portrait is organized by critical points (xq, yo) (i-e singular 
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Figure 3: Phase portrait of FRW model with dust and cosmo logical constant. 

solution for which yo = and dV/dx\(^xo,yo) = 0) and joining them trajecto- 
ries. Hence the phase portrait provides a global qualitative picture for the 
dynamics. Two phase portraits are equivalent if there exists an orientation 
preserving homeomorphism mapping integral curves of the systems. In our 
case the critical point for /3 + (1 + 3a) < is representing by static universe. 



where /3/(l + 3a) < 0. While constructing the phase portrait it is important 
for the system to be linearized near the critical points because the Hartman- 
Grobman theorem says that the original system is equivalent to its linear part 
in the nearby hyperbolic critical points (a critical point is hyperbolic (non- 
degenerate)) if there exists such i that ReAj 7^ 0, where Aj are eigenvalues 
of a linearization matrix. A character of critical points is determined by 
characteristic equation for the linearization matrix 



Since in our case TiA = only saddle points {d'^V/dx'^\(^xofl) < 0) or centres 
(c}^y^/c}x^|(a;o_o) > 0) are admissible. 

After some simple calculations one can check that a type of the critical 
point depends on the diagram of the potential, namely: 



Then 




(13) 




(14) 



detA 




3(l+a) 



(15) 
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Therefore for a > 1/3 there exists only one critical point which is of saddle 
type, i.e it is a structurally stable point. The full knowledge of the dynamical 
system comprises also its behaviour at infinity. To achieve this one usually 
transforms the phase space into a Poincare sphere. Then infinitely distant 
points of the phase plane are mapped onto the sphere's equator. The phase 
trajectories are mapped into corresponding curves on S"^. The character of 
critical points is preserved and new critical points representing asymptotic 
states of the system can appear at the equator. Then an orthogonal projec- 
tion of any hemisphere onto a tangent plane compactifies phase portrait. 



4 Horizon and flatness problem 

Due to existence of the first integral in such a system we can discuss some 
interesting properties of the evolutional path. There is a theorem ^] about 
nonexistence of particle horizon in such a model, namely if i? — and R < C, 
where C is a constant, then there is no particle horizon in the past. In terms 
of potentials, if V{R) —>■ constant for i? — > then there exists no particle 
horizon. It follows from the implication 

R<C >- — . (16) 

We can see from Fig. 1 and Fig. 2 in the S-J paper that all models with 
exception of case a) do not possess particle horizons. This is equivalent to 
nonexistence of the horizon problem in the corresponding models. However, 
it does not necessarily mean that if the condition R < C for i? — > is not 
satisfied then the horizon problem appears. In that case explicit evaluation 
of the integrals is required. We shall not deal with the problem in the present 
paper. 

As regards a flatness problem it does not appear in a cosmological model 
if for later times the curvature term does not dominate the matter term. In 
our model for the the problem not to appear we need the following condition 
to be satisfied 

W)l = IMt)l<^^. (17) 

Hence 

2 

Of course this is not the case considered either in the S-J or in the present 
paper, hence the flatness problem must be solved in a similar way as in the 
standard cosmology (e.g. using inflationary scenario). 
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5 Acceleration 



The case (3 < is particularly interesting because the corresponding cos- 
mological model without matter [A = 0) evolves with constant acceleration. 
With matter it accelerates always for later times independently of a. It 
follows from the relation 

r)V 1 A 2 

dR 2^ 2 ^ ' ^ ' 

Let us now differentiate both sides of the Eq. (4) with respect to r 

x = -- 5^(1 + 3ai)fiiox-(='"'+'). (19) 

If we want the model to accelerate at the present epoch we must put x > 
at X = 1. Hence we get the expected result 

^nonO > ^aO, (20) 

since this is the nonhomogeneity, which drives the acceleration. 



6 Magnitude-redshift relation 

It is well known that luminosity of observed objects, depends sensitively on 
the spatial geometry (curvature) and dynamics of the Universe. Therefore a 
cosmic distance measure, e.g. luminosity distance, depends on the present 
densities of different components filling up the universe and their equations 
of state. For this reason, the magnitude-redshift relation is proposed as a 
potential test for cosmological models and play important role in determining 
cosmological parameters [TE] . 

Let us consider an observer located at r = and at the moment t = to 
receiving light emitted at t from the source of absolute luminosity L located 
at the radial distance r. Cosmological redshift z of the source is related to t 
and to by the formula 

l + z _ Rjto) 

V{r,t) R{t)' ^ ' 

where the function V{r, t) is 

V{r,t) = l + -f3R{ty. (22) 
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If the apparent luminosity of the source measured by the observer is I, then 
the luminosity distance cIl of the source, defined 



is given by: 



dr 



(23) 



(24) 



Vir,t) • 

For historical reasons, the observed and absolute luminosities are defined 
respectively in terms of /^-corrected observed and absolute magnitudes m and 
M (/ = 10-2™/5 X 2.52 X 10-5ergcm-2s-2, L = IQ-^^/^ x 3.02 x lO^^ergs'^) 
|17j . When written in terms of m and M, Eq.(23) yields 



where 
and 



m{z,M,nao) = M + 51ogio[2^L(-2,f^ao)], 
M = M -5 logio Ho + 25 



(25) 

(26) 
(27) 



is a dimensionless luminosity distance in Mpc. 

The radial coordinate r in the expression (22) for di can be found by 
evaluation of the integral 

(28) 



X 

where x = R/Rq. Since in order to compare with observations the apparent 
luminosity m must be a function of the redshift z, the x parameter should 
be calculated from the relation 



z{x) = IH - 

X 4 



1 



dy 



^VL^qV^-'^^ + (1 - fi«o)l/^ 



(29) 



which follows from (21). Of course in a general case to find x analytically 
is not an easy task. Formally we have to find an inverse function x{z). The 
problem of finding dL{z) must be solved numerically and the simplest way to 
do that is to regard the function dL{z) as given in a parametric representation 



diix) 



U[z{x) + 1 - 1/a;] 



xHq 



(30) 
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and z{x) is given by (29). 

Similarly as in the paper of Jakacka and Stelmach we consider physically 
the simplest case, i.e. we assume that in the neighbourhood of the symmetry 
center the matter filling up the universe is a dust, hence a = 0. We define 
Q^o = Qao\a=o- In this case the expression (29) takes a form 



z{x) = 1 H - 

X 4 



dy 



^JVLm^y + (1 - ^mo)y^ 



(31) 



(^nono = 1 — ^mo), which cau be evaluated using Weierstrass elliptic function 
V pL4j- The invariants g2 and determining V function are 



92 = -^^moi^mO " 1), 93 = 0. 

Accuracy of the fit is characterized by the parameter 



X 



obs 

'0,i 



m, 



theor \ 2 
0,j / 



2 2 



(32) 



where m'^l is the measured value, 



j^theor jg value Calculated in the model 



^ J is a measurement error of m, while ^ is a measurement 



described above, 

error of z following from the dispersion in peculiar velocities of galaxies. 

We assume that supernovae measurements come with uncorrelated Gaus- 
sian errors and in this case the likelihood functions C can be determined from 
chi-squared statistics £ oc exp (— x^/2) [HIIH]- 



7 Results of the statistical analysis with the 
Perlmutter sample. 

In the paper of Stelmach and Jakacka it was shown that spherically symmet- 
ric Stephani cosmological model satisfying natural assumptions concerning 
local equation of state for the matter may explain supernovae data of Perl- 
mutter in the context of the magnitude-redshift relation. In the present paper 
we discuss this problem in more detail and we find numerical values of these 
parameters such that observational data are best fitted to our theoretical 
m — z curve. We take the whole sample of Perlmutter (all 60 galaxies) into 
account not rejecting any. Separately we analyse the sample of 54 super- 
novae (Perlmutter sample C), where 4 outliers and 2 reddened supernovae 
were excluded, finding no significant differences. In Fig. 4 we present plots 
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Figure 4: The Redshift-magnitude relation for Stephani model (Perlmutter 
Sample A). 

of magnitude-redsliift relations for the Stephani model and the Perlmutter 
sample A. 

The lower line corresponds to standard Einstein-de Sitter model, the 
line in the middle to our spherically symmetric Stephani model, i.e. Ai = 
—3.37, Qmo = 0.40 (this line is inseparable from the ACDM (Perlmutter) 
model with A = 0.7 and Qmo = 0.3), the upper line to Stephani model with 
Q^no = 0.3. From Fig. 4 we can see that the Perlmutter model and our best 
fit model are indistinguishable on the basis of the present available data. 

In Table 1 we present results of analysis with Ai obtained as a best fit 
for "classical" Perlmutter model with flrno = 0.3, Q\o = 0.7 (top line for each 
samples) and with marginalization over Ai (second line for each samples). 
The best fit has been obtained using Bayes technique. For Ai = —3.39, we 
obtain value of flmo = 0.37 (x^ = 96.25). Another good fit can be obtained 
with marginalization over Ai. In that case we obtain Ai = —3.37, Qmo = 
0.40. We test our results for the sample of 54 supernovae (Perlmutter sample 
C) and we realize that no significant differences were obtained. Note that in 
our model the numerical value for Q^o is larger than in the standard approach 
with cosmological constant, where flmo ~ 0.3 (fimo = 0.29 for sample A and 
n^no = 0.28 for sample C 14J). 

However, knowledge of the best-fit values alone has not sufficient scien- 
tific relevance, if confidence levels for parameter intervals are not presented 
too. Therefore, we carry out the model parameters estimation using the 
minimization procedure, based on the likelihood method. We could observe 
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that results obtained by this two methods are almost identical. At the con- 
fidence level of 68.3 % we obtain limits for values of parameters Ai and 
Qrno separately for samples A and C (Table 2). For example for Perlmutter 
sample A (with marginalization over Ai) we obtain that Ai = —3.37 and 
rimo = O-^Qtofs while for M = -3.39 we have fi„o = O.SYtof^ 

Varying ^2^0 and Ai we can find best-fit confidence regions in the {flmo, Ai) 
plane for our Perlmutter supernovae sample. In Fig. 5 we plot two confi- 
dence regions corresponding to levels of 95.5 (outer line) and of 68.3 (inner 
line) respectively. Since A4 is related to the absolute luminosity M and the 
Hubble constant Ho (see relation (26)), knowing M we can estimate Hq. For 
example assuming that 0.25 < Qmo < 0.40 (what is generally accepted) and 
M = -19.32 we find that 61.9 km/sMpc < Hq < 69.2 km/sMpc. 
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Figure 5: Confidence levels on the plane {flmo,-^)- (Perlmutter Sample A). 
The figure shows of the preferred value of Qmo and Ai at the confidence of 
level of 68.3 % and 95.4 



70. 



In Fig. 6 we show the levels of constant on the plane {^Imo, At). In this 
procedure we find the minimal value of x^, i-e., we consider best-fit values. 
The figure shows the preferred values of flmo and Ai. 

For a deeper statistical analysis of the Stephani model in explaining the 
currently accelerating universe we consider ID plot of the density distribution 
of l^nono- From this analysis one can obtain the limits at the la or 2a level. 
Fig. 7 shows the density distribution for finono in the Stephani model. This 
distribution is obtained from the marginalization over Ai. We obtain as a 
best fit value i^nono = 0.61. One can conclude that at the confidence level of 
68.3 % finono > 0.51 and finono < 0.70 while at the confidence level of 95.4 % 
we obtain that f2nonoe(0.41, 0.77). Fig. 8 shows the 1 dimensional density 
distribution for A4. 
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Figure 6: Levels of constant on the plane (flrno,-M) The figure shows the 
preferred value of (fimO;-^)- 




Figure 7: The density distribution for f^nono in the Stephani model (Perl- 
mutter Sample A), finono > 0.41, finono < 0.77 at the confidence level of 
95.4 %. 



13 



-3. '80 -3.' 60 -3. '40 -3. '20 -3. '00 

M 

Figure 8: The density distribution for M. in the Stephani model (Perl- 
mutter Sample A). The confidence level of 95.4 % is marked in the figure 
(A1e(-3.48,-3.25)). 

8 Statistical analysis with the Knop and Tonry 
samples. 

Because the Perlmutter sample was completed four years ago, it would be 
interesting to use newer supernovae observations. Lately Knop et al. [Zj 
have reegzamined the Permutter sample with host-galaxy extinction correctly 
applied. They chose from the Perlmutter sample these supernovae which were 
the more securely spectrally identified as type la and have reasonable colour 
measurements. They also included eleven new high redshift supernovae and 
a well known sample with low redshift supernowae. 

We have also decided to test our model using this new sample of super- 
novae. The mentioned authors distinguished few subsets of supernovae from 
this sample. We consider two of them. The first is a subset of 58 supernovae 
with extinction correction (Knop subsample 6; hereafter K6) and the second 
one a sample of 54 supernovae with low extinction (Knop subsample 3; here- 
after K3). Sample C and K3 are similarly constructed because both contain 
only low extinction supernovae. 

Another sample was presented by Tonry et al. [Sj who collected a large 
number of supernovae published by different authors and added eight new 
high redshift SN la. This sample of 230 Sne la was recalibrated with con- 
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sistent zero point. Whenever it was possible the extinctions estimates and 
distance fitting were recomputed. However, none of the methods was able 
to apply to all supernovae (for details see Table 8 in jH]). This sample was 
improved by Barris who added 23 high redshift supernovae including 15 at 
z > 0.7 doubling the published number of object at this redshifts [6J. 

Despite of the mentioned above problems, the analysis of our model using 
this sample of supernovae could be interesting. We decide to analyse four 
subsamples. First, the full Tonry/Barris sample of 253 SNe la (hereafter 
sample TBa) is considered. The sample of 218 SNe la (hereafter sample 
TBb) consists of low extinction supernovae only (median V band extinction 
Ay < 0.5). Because the Tonry sample has a lot of outliers especially in low 
redshift, we separately analysed the sample where all low redshift {z < 0.01) 
supernovae are excluded. This sample again contains 218 SN la, but they 
are different than that belonging to the sample TBb (hereafter sample TBc). 
In the sample of 193 SN la all supernovae with low redshift and and high 
extinction are omitted (hereafter sample TBd). 
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Figure 9: The density distribution for i^nono in the Stephani model 
(Tonry/Barris sample TBa). f2nonoe(0.52, 0.68) at the confidence level of 
95.4 %. 

Tonry and Barris [3 IE] presented redshift and luminosity distance obser- 
vations for their sample of supernovae. Therefore, Eqs. (25) and (26) should 
be modified [T^ : 

m - M = 5 logio(I^L)Tonry " 5 log^Q 65 + 25 (33) 
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and 

= -51ogioi^o + 25. (34) 

For Ho = 65 km Mpc"^ we obtain M = 15.935. 

The results obtained with the new sample is very similar to that obtained 
with Perlmutter sample but now errors decrease. The results with Tonry 
sample are almost identical to that with Perlmutter sample. For the sample 
TBa we obtain that Qmo = 0.39, while for ACDM model flmo = 0.32 (for 
the sample TBd ^2^0 = 0.33 0). Using Tonry sample we obtain at the 
confidence level of 95 % hmit for nnone(0.52, 0.68) (see Fig. 9). However with 
the Knop sample we obtain the value flmo niore closed to value 0.3, what was 
obtained from CMBR and extragalactic data |2m 12 Ij. than in the previous 
case. However, it should be noted, that for ACDM model we obtain for the 
sample K3 Q^o = 0.25 while for the sample K6 Qmo = 0.28. It means, that 
our conclusion that the numerical value for Qmo in our model is larger than 
in the standard approach (ACDM model) is still valid. 
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Figure 10: The Redshift-magnitude relation for the Stephani model (Knop 
Sample K3). 

In Fig. 10 we present plots of redshift-magnitude relations for Stephani 
model Perlmutter sample A. The lower line corresponds to the standard 
Einstein-de Sitter model, the line in the middle to our spherically symmetric 
Stephani model, i.e. j\4 = —3.46, Qmo = 0-32 (this line is inseparable from 
the ACDM model with A = 0.75 and Qmo = 0.25), the upper line to Stephani 
model with Qmo = 0.25. From Fig. 10 we can see that the ACDM model and 
our best fit model are indistinguishable also on the basis of the Knop sample. 



16 




0.00 0.20 0.40 0.60 0.80 1.00 1.20 



Figure 11: The density distribution for r2nono in the Stephani model (Knop 
sample K3). f2nonoe(0.55, 0.79) at the confidence level of 95.4 %. 
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Figure 12: The density distribution for M. in the Stephani model (Knop 
Sample K3). Ale(-3.54, -3.37)) at the confidence level of 95.4 %. 
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Using the Knop sample (Figs 11 and 12) we obtain finone(0.55, 0.79) at 
the confidence level of 95 %. Because Knop discuss very carefully extinction 
correction and as a result his sample has extinction correctly applied, we 
think that using the limit obtained from the Knop's sample is the most 
appropriate. Our results show that Stephani model is consistent with SNIa 
data at the 95 % confidence level. Our study shows that Stephani model is 
very good fit to latest supernovae data and should be treated as a alternative 
to ACDM (Perlmutter) model. 



9 CMB peaks in the Stephani model 

The CMB peaks arise from acoustic oscillations of the primeval plasma. 
Physically these oscillations represent hot and cold spots. Thus, the wave- 
length of the perturbation which contributes the most to the density distri- 
bution at the time of the last scattering corresponds to a peak in the power 
spectrum. In the Legendre multipole space this corresponds to the angle sub- 
tended by the sound horizon at the last scattering epoch. Higher harmonics 
of the principal oscillations, which oscillated more than once, correspond to 
secondary peaks. 

It is well known that the locations of the peaks are very sensitive to the 
variations in the parameters of the model. Therefore, it can be used as a 
sensitive probe to constraint the cosmological parameters and discriminate 
among various models. 

The locations of the peaks are set by the acoustic scale Ia which can be 
defined in terms of an angle 6 a subtended by the sound horizon at the last 
scattering surface. In a general case calculation of the angle 6 corresponding 
to acoustic or particle horizon existing at the recombination epoch is not an 
easy task even in the Friedmann model. One of the way to find 6 is to solve 
an appropriate Euler-Lagrange problem (|221, |IZI)- If the angle 9a is small 
(in general it is small, because Cg « c), hence the acoustic scale Ia = t^/Oa 
is given by 

jj^ dx/h{x) 
Jq Cs[x)dx I h[x) 

where 

h{x) = a/ ilro + ^rnOX + i^nonOa:^ (36) 

and the relation between z and x is given by the equation (31). Cs{x) is 
a speed of sound in the plasma and varies with the expansion (we assume 
additionally presence of radiation in the model and Qro is an energy density 
parameter corresponding to radiation at the present epoch). Similarly flmo 
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corresponds to nonrelativistic matter and finono to nonhomogeneity. The 
sound velocity can be calculated from the formula 



= ' , oo ^ o . (37) 



In the case without cosmological term the parameters Q^no and finono are not 
independent, and finono can be expressed as 

^nonO = 1 ~ ^rO ^ ^mO- (38) 

In the model of primeval plasma, there is a simple relation 

Im ~ him - (j)m) (39) 

between the location of the m-th peak and the acoustic scale I24 j . The 
prior assumptions in our calculations are flro = 9.89 ■ 10~^, flbo = 0.05 
(baryonic matter energy density), and the spectral index for initial density 
perturbations is n = 1. Moreover we put for the present value of the Hubble 
parameter Hq = 65 km/sMpc. 

The phase shift is caused by the pre-recombination physics (plasma driv- 
ing effect) and hence, is not significantly influenced by the Stephani term 
at that epoch. Therefore the phase shift (pm can be taken from standard 
cosmology 



0.267 



0.3 



0.1 



(40) 



where Qboh"^ = 0.02, r(zdec) = Pr{zdec) / Pm{zdec) = ^ro/^moXdec- Radiatiou 
energy is composed of two components: electromagnetic energy and neutrino 
energy flro = ^-yo + ^uo, ^70 = 2.48/i"^ ■ 10"^, Q^o = ■ 10"^. r{zdec) is 

the ratio of radiation to matter densities at the surface of last scattering. 

In the Friedman models (finono = 0) with possible presence of cosmological 
constant we obtain: 

^mO = 1-0 f^nonO = 0.: ^peak.l = 203, /peak,2 = 471, /peak,3 = 739, 
^mO = 0.3 = 0.7: ^peak.l = 225, /peak,2 = 536, /peak,3 = 847. 

The influence of f^nono on the location of the peaks results in shifting them 
towards higher values of / in comparison to the model with Qro- For example, 
for flbo = 0.05, h = 0.65, the different choices of finono yield 

^mO = 0.4 finonO = 0.6: /poak,l = 212, /peak,2 = 500, /peak,3 = 789, 
^mO = 0.3 finonO = 0.7: /peak,l = 216, /peak,2 = 514, /peak,3 = 812. 

On the other hand from the Boomerang observations [21] we obtain 
^peak,i = 200 -7- 223, /peak,2 = 509 561. We also compare the results 
from our model to recent bounds on the location of the first two peaks 
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obtained by WMAP experiment 123 • Namely /peak.i = 220.1 ± 0.8, 
^peak,2 = 546 ± 10, together with the bound of the location of the third 
peak obtained by Boomerang experiment /peak,3 = 825ti3 which lead to quite 
strong constraints on the model parameters. These constraints can be sum- 
marized as follows. The Stephani model is in agreement with observations 
and we conclude that the influence of the term f^nono is not very significant 
in our case. However, phase shift is taken from the standard cosmology, 
i.e., we assume that the contribution from the "Stephani nonhomogeneity" is 
insignificant at the pre-recombination epoch. If this assumption is not valid 
then the limit from CMB will change. 

It should be also noted that the Big-Bang Nucleosynthesis (BBN) is a very 
well tested area of cosmology and does not allow for significant deviation 
from the standard expansion law, apart from very early times before the 
onset of the BBN. The consistency with the BBN seems to be necessary 
in the Stephani model. We consider only non-relativistic matter p in the 
Stephani model. In this case we could approximate that f^nono term scales 
like (1-1-2;). It is clear that the contribution of the "Stephani nonhomogeneity 
term" cannot dominate over the standard radiation term before the onset of 
BBN, i.e., for z = 10^ and as a result has no effect on the BBN. 

10 Conclusions 

Our investigations show that the Stephani model is an excellent fit to both 
Perlmutter data points and currently available Knop's data points. Sum- 
marizing, relatively large value of finono = O.GI^qjq can explain supernovae 
data of Perlmutter. In the range of observed redshifts of supernovae {z < 1) 
the curve corresponding to m — z relation in our model is almost indistin- 
guishable from the ACDM (Perlmutter) model with cosmo logical constant. 
With future data from SNAP the error bars in the extimation of the model 
parameters will be reduced significantly. Together with new limits for Qmfl 
(obtained from extragalactic data) it will be possible discriminate between 
Stephani model and ACDM (Perlmutter) model. 

Locally our spherically symmetric Stephani universe is indistinguishable 
from the Friedman one (the same equation of state, the same local geome- 
try), hence local physical cosmology in both models must be the same. In 
other words standard scenario of the evolution of the universe remains valid. 
Especially in the early universe effects of non-homogeneity become negligible. 
This is easily seen if we evaluate the energy density parameter for matter 
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ClS db function of the scale factor R 



1 



(41) 



1 - 3/3C^i?2+3a ■ 



In the early universe relativistic matter dominates {a = 1/3) and from the 
above relation it follows that ^^(-R) tends very quickly to 1 for R ^ 0. It 
means that flrion{R) — > 0. This justifies application of standard physical cos- 
mology (including e.g. primordial nucleosynthesis) for our nonhomogeneous 
model. 

Freese and Lewis in I^Hl have recently proposed interesting alternative 
model explaining the currently accelerating Universe. In this model which 
is called the Cardassian model the standard Friedmann-Robertson- Walker 
(FRW) equation is modified by the presence of an additional term p", namely 



where p is a total energy density (matter + radiation), and B is a positive 
constant. This proposal seems to be attractive because the expansion of the 
universe is accelerated automatically due to the presence of the additional 
term (if we put B = then the standard FRW equation is recovered) without 
postulating existence of unknown form of dark energy. 

Note that putting p = 3pmo/R^ (and hence a = 0) and n = 1/3 and 
B = —3P/^y3pmo the above equation could be rewritten in the form: 



Comparing to (1) we realize that we obtain correspondence between Cardas- 
sian and Stephani models. 

When for simplicity we assume that the energy density parameter for 
radiation matter vanishes {Qro = 0) we recover the model analysed in our 
previous paper [HO]- We note that considered by us spherically symmetric 
Stephani model is a special realization or in other words is dynamically equiv- 
alent to the Cardassian model for n=l/3. It is interesting that this value of 
n parameter is reasonable from the observational joint test analysis of both 
CMB from and SNIa 

Let us note that in our model there is the cosmic coincidence problem, i.e. 
problem why did the nonhomogeneity start to dominate the present evolution 
of the Universe only fairly recently? There is no satisfactory solution of this 
problem and some other face of fine tuning coincidence is required. 



H' = f^ + Bp-, 



(42) 




(43) 



21 



References 

Jakacka I and Stelmach J 2001 Class. Quantum Grav. 18 2643 
Krasinski A 1985 GRG 15 673 

Krasinski A 1997 Inhomogeneous Gosmological Models Cambridge Uni- 
versity Press Cambridge 1997 p. 317 

Perlmutter S et al. 1999 Astrophys. J. 517 565 

Tonry J L et al. 2003 Astrophys. J. 594 1, jastro-ph/0305008l 



Harris B J et al. 2003|astro-ph/0310843 



Knop RAet al. 2003, |astro-ph/ 0309368 

Perlmutter S et al. 1998 Nature 391 51, Perlmutter S et al. 1999, Aas- 
trophys. J. 517 565, Garnavich P M et al. 1998 Astrophys. J. Lett. 493 
L53, Riess A G et al. 1998 Astron. J. 116 1009 

Caldwell R R, Dave R and Steinhardt P J, 1998 Phys. Rev. Lett. 80 
1582, Zlatev I, Wang L and Steinhardt P J, 1999 Phys. Rev. Lett. 82 
896, Steinhardt P J, Wang L and Zlatev I 1999 Phys. Rev. D59 123504 

Dabrowski M P and Stelmach J, in Large Scale Structure in the Universe 
lAU Symposium No 130 eds. Jean Adouze, Marie-Christine Palletan and 
Alex Szalay, (Kluwer Academic Publishers, 1988), 566. 

Dabrowski M P and Stelmach 1987 J Astron. Journ. 93 1373 

Stelmach J, D§browski M P and Byrka R 1993 Nuclear Physics B406 
471 

Celerier M N 2003 Astronomy and Astrophysics 353 63 

D§browski M P and Stelmach J 1986 Ann. Phys (N. Y.) 166 422 

Szydlowski M, Krawiec A and Dabrowski M P 2002 Phys. Rev D66 
064003 

Dabrowski M P and Hendry M A 1998 Astrophys. J. 498 67 
Weinberg S 1972 Gravitation and Gosmology (Wiley, New York, 1972) 
Riess A G 1998 Astron. J. 116 1009 



22 



Williams B F et al. 2003|astro-ph/0310432 



Peebles P J E, Ratra B 2003 Rev. Mod. Phys. 75 559, |astro-ph/0207347 
Lahav O 2002 |^itro-pli/0208297| 

Narlikar J V 1983 Introduction to Cosmology, (Jones and Barlett, Boston 
1983) 

Doran M, Lilley M, Schwindt J and Wetterich C 2001 Astrophys. J. 559 
501 

Hu W, Fukugita M, Zaldarriaga M and Tegmark M 2001 Astrophys. J. 
549 669 

Ichiki K, Yahiro M, Kajino T, Orito M and Mathews G J 2002 Phys. 
Rev. D 66 043521 

de Bernardis P et al. 2002 Astrophys. J. 564 559 
Spergel D N et al. 2003, |astro-ph/0302209| 



Page L et al. 2003, |astro-ph/ 0302220 

Freese K and Lewis M 2002 Phys. Lett. B 5401, astro- ph/020122 9| 

Godlowski W, Szydlowski M and Krawiec A 2004 accepted for publica- 
tion in Aastrophys. J., astro-ph/0309569, 

[31] Sen S 2003 AA588 1 |astro-ph/0211634| 



23 



Table 1: Results of the statistical analysis of the model using the best fit 
with minimum (denoted by BF). The case in which we marginalize over 
M. is denoted by /A. 



sample 




^"^iioiiO 


M 




method 


A 


0.37 


0.63 


-3.39 


96.3 




RP 

JJ J. 




0.40 


0.60 


-3.37 


96.1 


A/f 


RF 


C 


0.36 


0.64 


-3.42 


53.4 




RF 




0.36 


0.64 


-3.42 


53.4 


A/f 


RF 


K3 


0.29 


0.71 


-3.48 


61.3 




BF 




0.32 


0.68 


-3.46 


61.0 


M, 


BF 


K6 


0.36 


0.64 


-3.53 


56.3 




BF 




0.36 


0.64 


-3.51 


56.0 


M, 


BF 


TBa 


0.38 


0.62 


15.905 


262.3 




BF 




0.39 


0.61 


15.915 


262.2 


M, 


BF 


TBb 


0.39 


0.61 


15.925 


204.8 




BF 




0.39 


0.61 


15.925 


204.8 


M, 


BF 


TBc 


0.39 


0.61 


15.915 


229.3 




BF 




0.41 


0.59 


15.925 


229.2 


M, 


BF 


TBd 


0.39 


0.61 


15.925 


192.7 




BF 




0.39 


0.61 


15.925 


192.7 


M, 


BF 
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Table 2: Results of the statistical analysis of the model using the likelihood 
method (denoted by L). The case in which we marginalize over M. is denoted 
hyM. 



sample 






method 


A 


n 07+U.U5 
'-'■"5' -0.05 


-3.39 


L 




u.oy_o,o8 


_o oq+0.06 


LM 


C 


'-'•'J"-0.05 


-3.42 


L 




'-'•'^"-0.08 


q /1 9-1-0.05 
"-'•^^-0.05 


LM 


K3 


"J-^^-o.os 


-3.48 


L 




'-'•^^-0.06 


_o 4f5-H0.05 


LM 


K6 


'-'■'^-'--0.06 


-3.53 


L 




n qq+0.08 
'-'•'-''-'-0.07 


q p: 1-^0.03 
'J -^J^ -0.03 


LAI 


TBa 


0.38l°:°i 


15.905 


L 




0.39ini 


15.915l°:2^ 


LAI 


TBb 


0.39lH^ 


15.925 


L 




0.39iHl 


15.925^Q 02 


LAI 


TBc 


n Qq+0.02 

'-'•'-'^-0.02 


15.915 


L 




'-'•^U-0.04 


15.925^Q 02 


LAI 


TBd 


0.391°:°^ 


15.925 


L 




0.391^^ 


15.935^0 02 


LM 
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